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ABSTRACT 

We present the results of a kinematical study of the W 50 nebula using high resolution radio 
observations from the Very Large Array (VLA) spanning a 12-year period, sampled in 1984, 
1993 and 1996. We conduct a careful analysis of the proper motions of the radio filaments 
across the W 50 nebula at each epoch, and detect no significant motion for them during this 
period. The apparent lack of movement in the radio filaments mandates either (i) a high de- 
gree of deceleration of SS 433's jet ejecta in the W 50 nebula, or (ii) that the lobes of W 50 
formed a long time ago in SS 433's history, during a jet outburst with appreciably different 
characteristics to the well-known precessing jet state observed in SS 433 at the present day. 
We discuss the possible scenarios which could explain this result, with relevance to the nature 
of SS 433's current jet activity. 

Key words: stars: individual:SS 433, ISM: kinematics and dynamics, ISM: jets and outflows, 
ISM: supernova remnants, hydrodynamics 



1 INTRODUCTION 

The spectacular interaction between the relativistic jets of the mi- 
croquasar SS433 and the supernova remnant (SNR) W 50, has been 
captured at several epochs by radio telescopes that are sensitive to 
both its angular extent (approximately two degrees across the sky) 
as well as to its fine structure on arcsecond scales. The connection 
between microquasar and SNR has been confirmed, as the distance 
to this structure is now very well known via two distinct methods. 
In the first instance, a distance to SS433 of 5.5 + 0.2 kpc was de- 
termined by examination of light-travel time induced aberrations 
of the arcsec-scale jets of SS 433 ( [Blundell & Bowler|2004] ). The 
second method is based upon the detection of neutral hydrogen, in 
both absorption and in emission, at a high Galactic rotation speed 
corresponding to the same distance of 5.5 kpc to W50 (Lockman 
|et al.||2007| ). Accurate distances to objects within our Galaxy are 
still rare, but when such objects are identified their angular sizes 
translate directly into physical sizes, proper motions convert nicely 
to velocities, and this information enables a quantitative assessment 
of the system energetics. 

The central engine in SS 433 is known to eject jet material 
with a mean speed of ~0.26c as determined from optical spec- 
troscopy (Eikenberry et al. 2001) and this is found to be consis- 
tent with proper motions of the jet knots using long-baseline radio 
interferometry (Vermeulen et al. 1987, Paragi et al. 19 99} with a 
standard deviation in the jet speed of 0.01c ( [Blundell & Bowler| 
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2004). The recent review of|Pabrika ( 2004) provides a detailed de- 
scription of SS 433's behaviour and a comprehensive summary of 
different types of observation of this system. 

Although this paper focuses upon radio observations, evidence 
for the SS 433-W 50 interaction has been reported in several other 
wavelength regimes. Optical maps of W50 ( [Zealey et al."l[T980l 
|Boumis et al.| [2007) show bright filaments approximately 30' to 
the east and west of SS 433 that are suggested to be excited by the 
relativistic beams of SS433. Broadband X-ray observations ( |Safi-| 
Harb & Oegelman 1997; Brinkmann et al. 2007) show emission 
peaking just radially inwards of the extremities of the radio protru- 
sions (e.g. |Blundell & Hirst 2011 fig. 3) with hard non-thermal 
X-ray emission at 35' from SS433, and a general softening of the 
X-ray spectrum with increasing distance from SS433. The X-ray 
emission originates near to the mean precession axis, well within 
the 20° semi-angle cone traced out by the precessing jets. This has 
led to suggestions that SS 433's jets are being refocused back to- 
ward the mean jet axis ( [Peter & Ei chler 1993), or that the X-ray 
emission arises due to reflective hydrodynamic shocks propagat- 
ing along the symmetry axis because of the SNR-jet interaction 
( Velazquez & Raga 2000). However, more recent hydrodynamic 
simulations of the SS 433-W 50 interaction show a new mechanism 
for hydrodynamic refocusing of conical jets (Goodall et al.|20li] ), 
where W50's nebular morphology and SS 433's refocused jets are 
both reproduced to a good approximation. X-rays have been de- 
tected from SS 433 itself by the Chandra satellite; on the basis of 
some of these observations Migliari et al. ( 2002 ]) suggest in- situ re- 
heating of the baryonic components, while Migliari et al. ( 2005) 
interpret the variability of the arcsec-scale X-ray jets as evidence 
of a second, faster outflow connected with the jets. Variability and 
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deviation from the standard kinematic model for the radio and op- 
tical jets in SS 433 have been analysed and considered by Blundell 
|& Bowler| ( [2005l > and |Blundell eTall { [20071 ). 

Whilst the kinematics of SS433's jets have been generally 
well- studied near the point of jet-launch at both optical and radio 
wavelengths, the relativistic propagation of jet material at the outer- 
most periphery has not yet been observed. Indeed, any deceleration 
of the jet material, or the point at which deceleration happens, has 
yet to be established. 

The SS433-W50 system is amenable to numerical simula- 
tions, and the results of a wide range of hydrodynamic models 
are presented in a companion paper (Goodall et al. 2011). In the 
present paper, we present a detailed examination of three VLA ra- 
dio images of W 50 spanning 12 years. We aim to ascertain whether 
meaningful constraints can be placed on the proper motion of the 
persistent radio features (the filamentary structure in W 50) where 
the jets have impinged and penetrated through the supernova rem- 
nant (SNR) shell. In particular, we aim to identify any evidence for 
the deceleration of the jets at the furthest point they are observed 
from the microquasar nucleus. 



2 DATA ANALYSIS 
2.1 VLA Data Reduction 



For the datasets used in this analysis, the image published by |Dub-| 
|ner et al.| fl998) was kindly made available to us by Gloria Dubner, 
and the NVSS image was downloaded from the NVSS "postage 
st amp server" ([Condon et al.|1998| >. The data originally observed 
by |Elston & Baum| ( |1987| ) were retrieved from the NRAO VLA data 
archive. These data comprised nine telescope pointings, observed 
at 1.4 GHz in the VLA's most compact configuration, D-array, and 
were reduced individually using standard procedures in AIPS. The 
individual pointings were then sequenced together using bespoke 
code written in PerlDL. For clarity we refer to these observations 
by source and observation year according to their respective FITS 
headers, as Elston-84 (date-obs: 12/08/1984), NVSS-93 (date-obs: 
15/11/1993), and Dubner-96 (date-obs: 19/08/1996) respectively. 



2.2 Image Preparation 

Three epochs of W50 observations from the VLA in L-band and 
D-array configuration, provide us with maps of very similar resolu- 
tion (the VLA synthesised beamwidtlQin this configuration is 44") 
and size, with differences being due to the restoring beam used for 
each observation and the fact that each image is comprised of a se- 
ries of mosaicked pointings. An interactive kinematics p ipeline was 
written by PTG using the Perl Data Languag^](PDL, Glazebrook 
& Economou (1997)) to extract the accurate fitted coordinates of 
the SNR filaments from the radio observations. The pipeline per- 
forms a series of image preparation steps to ensure accurate image 
matching before kinematics can be measured, and the basic stages 
are outlined below. 

The images were scaled-up from their original sizes (by a fac- 
tor of at least 5) such that the image sizes and pixel scales ap- 
proximately matched, using a high-order polynomial interpolation 



1 http://www.vla.nrao.edu/astro/guides/vlas/current/node 1 0.html 

2 The Perl Data Language (PDL) has been developed by K. Glazebrook, 
J. Brinchmann, J. Cerney, C. DeForest, D. Hunt, T. Jenness, T. Luka, R. 
Schwebel, and C. Soeller and can be obtained from http://pdl.perl.org. 



method to ensure accurate resampling. Each image was then nor- 
malised such that the brightest pixel (always centred upon the unre- 
solved source SS 433) has a value of unity. An accurate astrometric 
correspondence between the three observations was created using a 
select group of 15 bright point sources common to all three images 
(see Figures B.l to B.3 in Appendix B). This was done by fitting 
Gaussians to the point source centroids, and each fit was inspected 
by eye to ensure positional reliability. The averaged coordinates 
C=(x c ,y c ) of all 15 point sources (the centre of mass if all point 
sources are assigned an equal mass) was chosen as the principal 
reference point for each image: 
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where n is the number of point sources used, and the error asso- 
ciated with any one measured centroid position is therefore min- 
imised. The 6-character capitalised keywords (CRPIX1 etc) refer to 
the standard FITS header keywords defined by W ells et aL] ( |1981| ). 
We define the east centre E-(x E , y E ) as the average coordinates of 
all the point sources to the east of SS 433 (point sources 2-9), and 
west centre W=(x w , j w ) as the average coordinates of the remain- 
ing point sources including SS 433 (point sources 10-15 and 1) such 
that: 
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where a = 8 and b = 7 are the number of sources used for the east 
and west sides respectively. The separation between the east centre 
(E) and west centre (W) then defines an accurate image pixel scale 
for each image: 
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The angle 6 made by the vector eW (ioinins the east and west 
centres) with the horizontal axis is indicative of any residual rota- 
tion that must be performed to match the images: 
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The Elston-84 and NVSS-93 images were scaled and rotated 
to match the new header values of the Dubner-96 image, and the 
Gaussian fitting to the 15 point source centroids was repeated to 
establish the correct header information for the newly transformed 
images, followed by a convolution to match both of Elston-84 
and NVSS-93 to the slightly coarser resolution of # res =56"of the 
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Figure 1. An example of the Gaussian centroid fitting to determine 
the filament mean positions along its length. The filament has been 
rotated such that its longest axis is along the vertical, and the Gaus- 
sian curves are then fit horizontally line by line, before rotating 
the coordinates back to the image frame. Note that only every 25th 
sampling is shown, and only every second datapoint per intensity 
profile. Difficulties arise when the structure of the filament is more 
extended, causing the horizontal intensity profile to diverge from a 
Gaussian shape. To compensate for this, each Gaussian curve is fit 
individually to determine a reliable range for the fit, and each fit is 
inspected by-eye to ensure a good match to the peak-position (the 
wings of the Gaussian are not important here). The white points 
represent the data, and the grey lines indicate the best fit to each 
line of data shown. Each Gaussian curve here has been overlaid 
upon the image, so that its peak aligns with the y-position from 
which it was taken. 



Dubner-96 image. Finally, the images were matched to within 0.01 
pixels by performing an image product maximisation based upon 
the 15 point sources. A final transformation was applied using a 
high-order polynomial interpolation, and the image FITS headers 
were updated one last time by fitting Gaussian centroids to the point 
sources in the matched image^] 

Throughout the analysis, non-linear corrections in astrometry 
were not accounted for, because the linear approximation is satis- 
factory over the degree scales involved. 



2.3 Velocity resolution 

The bright filaments in the images generally have a relatively high 
signal-to-noise and the Gaussian curves we fit to the cross sections 
of the filaments are well sampled across a given reliable range 
(determined by inspection, detailed in j j2.4| l. Under these condi- 
tions, the accuracy with which one can determine the centroid po- 
sition from a Gaussian fit is typically much finer than the Gaussian 
FWHM, often by a factor of 10 or more depending on the signal- 
to-noise. Hence, the Gaussian centroid-fitting accuracy x res can be 
written as 



3600 180 



</ss433 = 0.149 pc = 0.487 ly 



(11) 



where dss433 = 5500 pc is the distance from us to SS 433 and W 50, 
and //fit = 0.1 represents the modest factor of 10 times better ac- 
curacy in determining the centroid position, as compared with the 
image resolution. The time baseline £ b ase between the Elston-84 and 
Dubner-96 images is 12 years and 4 days (with three of those years 
being leap years) giving a total of 4387 days, or of ? b ase = 12.01 
years based upon the 3 65. 25 -day year. The velocity resolution in 
units of the speed of light is then: 



0.0405 c. 



(12) 



This velocity resolution is a factor of 6.5 times smaller than 
the measured speed of SS433's jets of Vj et = 0.2647c (Eikenberry 
|et al.|200l) close to their launch point, and so the kinematics of 
the movement of filaments associated with the jets (or turbulence 
caused by the jets) should be measurable, unless significant decel- 
eration has occurred. 



2.4 Measuring Filament Kinematics in W 50 

The ten brightest filament-like structures common to all three im- 
ages were selected and labelled A-H as shown in Figures |B1[ |B2[ 
and |B3| Although other common filamentary features are present 
between the three images, only the ten brightest were used based 
upon the criteria of having a a signal-to-noise ratio greater than five, 
and being at least four times as long as wide, in all three images. 
Also, the three images are not of the same quality; the Dubner-96 
image is by far the deepest, followed by the Elston-84 image. NVSS 
was a survey designed to map the bright radio sky quickly, and so 
the NVSS-93 image of W 50 is much lower fidelity due to the lower 
integration time on target. For this reason, the filament brightnesses 
vary considerably across the three images, and the NVSS-93 image 



3 An animation sequencing the final images is available here: 
http : / lwww-astro.physics.ox.ac.ukl~ptglRESEARCHIresearch.html 



© 0000 RAS, MNRAS 000, 000-000 



4 R T. Goodall, K. M. Blundell and S. J. Bell Burnell 



Displacement from SS433 (arcminutes) 
-70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 




I i I i I i I i I i I i l i I i I i I i l_ 

-100 -80 -60 -40 -20 20 40 60 80 100 

Displacement from SS433 (pc) 




Figure 2. (a): The filament positions of each of the three observations are overlaid upon a contour map of the Dubner-96 observation. The 
yellow-black dashed lines represent the pseudo-dataset generated by projecting the Elston-84 filament positions 12 years into the future, 
assuming a speed of 0.26c in the direction radially away from SS 433. The brown arrow indicates SS 433 's jet precession axis, which is also 
W50's symmetry axis. The magenta arrows in dicate the edges of SS 433's precession cone in the plane of the sky, based upon the cone 
semi-angle at the current epoch of 0j et = 20.92° (Eik enberry et al.|2001| ). The points P A , Pb, Pc. and P D indicate the locations where SS 433's 
current jet configuration would penetrate the W 50 nebula. Note the remarkable alignment between the bright filament nearest to SS 433 and 
the north-east jet cone line at image coordinates ~(-10,+3) parsecs. The radial distance d A from SS433 to point P A is approximately 50 
parsecs. The errors in the fitted positions of the filaments are of the order of (or smaller than) the thickness of the filament lines plotted on 
the figure, (b & c): A zoom-in on each of the eastern and western filaments. Note that filaments A and J which are the most prominent from 
the radio images, are almost perpendicular to the precession axis. 



was the limiting factor in choosing suitably bright filaments. Fil- 
ament B was later rejected from the NVSS-93 dataseQ during the 
analysis by the pipeline, as not being bright enough determine its 
position accurately. Using the interactive kinematics pipeline (as 

4 The lack of Filament B from the NVSS observation does not affect any 
of the conclusions of this paper. 



described in |2.2| ), the longest axis of each filament was determined 
and rotated to the vertical, such that the Gaussian curve fitting could 
be applied to the horizontal direction, to each pixel-row of the fil- 
ament region. The fitting process could not be automate^ due to 
extended structural variations in the background emission around 

5 This involved over 2000 Gaussian fit curves per image. 
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each filament, and each Gaussian curve fit was inspected by-eye 
to ensure the best-fit was measured for the peak position of the 
filament and not to background structure which contributed to the 
Gaussian wings. An example of the Gaussian fit curves to Filament- 
A is given in Figure [T] and the complete set of filaments from each 
image is given in appendix B. 



3 RESULTS 

The filament positions for each of the three epochs are overlaid 
together upon a contour image of the Dubner-96 observation in 
Figure [2] A pseudo-dataset was generated using the data from the 
Elston-84 observation, by projecting those filament positions 12 
years into the future (based upon SS 433's current jet speed at the 
point of launch) in order to compare the expected pseudo-filament 
positions with the Dubner-96 data. This projection was only ap- 
plied to filaments that fall within SS 433's precessional jet cone, as 
material far from the mean jet axis will not have interacted directly 
with the jets and cannot be expected to propagate at similar speeds. 
Inspection of Figure [2]re veals that the filament positions from each 
of the three observations are remarkably coincident, with little in- 
dication of any conclusive movement, whereas the projected fila- 
ments are noticeably displaced from the rest. This raises two possi- 
bilities: 

(i) either considerable deceleration of the jet ejecta has occurred 
between jet-launch and the filaments of the jet lobes, or 

(ii) the filaments in W50 have not interacted with the current 
jet (of speed ~0.26c) observed to be active near the central engine 
SS433. 

A velocity map (Fig. [3} taken from high resolution hydrodynamic 
simulations of the SS433-W50 interaction (detailed in Goodall 
|et al.| ( |2011| )) shows no signs of the required deceleration along 
the jet axis, and therefore supports the second possibility. 



4 DISCUSSION 

The analysis reveals no discernible changes in the filament posi- 
tions during the 12-year baseline covered by these observations. We 
therefore infer an upper limit to the proper motion of the filaments 
equal to the velocity resolution of this investigation: v res = 0.0405c 
or// = 0.466" yr -1 - 

One conclusion to draw from this (as from (i) above) is that 
considerable deceleration of the jet ejecta has occurred through in- 
teractions with the interstellar medium (ISM) present along the path 
of the jet lobes of W50, causing a reduction in the speed of the ma- 
terial in the lobes of W50. This deceleration must be sufficient to 
reduce the speed of the jet material to below the detectable veloc- 
ity resolution (v res = 0.0405c) of this kinematic study. These data 
would imply a lower limit of ~ 85% kinetic energy loss in the jets, 
and this energy would have been transferred to the ISM of the lobes. 
However, high resolution hydrodynamic simulations of the jets in 
SS433 (Goodall et al. 2011) reveal that, for any reasonable com- 
bination of jet mass-loss rate in the range 10~ 5 - 10~ 4 M yr -1 and 
ISM density in the range 0. 1 - 1 particles cm -3 , a persistently active 
jet of speed v jet = 0.26c does not suffer such extreme deceleration 
(see Figure [3}. We stress however, that this model is based upon a 
persistent jet, whereby kinetic energy is continually injected from 
the central engine in SS 433 to replenish the kinetic energy (and 



momentum) lost to the jet, and transferred to the surroundings via 
interactions with the ISM. 

In order to create a possible scenario that is in agreement both 
with the kinematics data presented here and the hydrodynamical 
simulations of the companion paper Goodall et al. ( 201 1 ), we must 
invoke the idea of intermittent, rather than persistent, jet behaviour 
in SS433. The case for discontinuous jet activity in SS 433's his- 
tory can be realised through assimilation of the relevant facts: 

(i) Our hydrodynamic simulations ( Goo dall et al.|201l) follow 
the evolution of a supernova explosion, until the SNR reaches 45 pc 
in radius, equal to that of the circular region of the W 50 nebula. The 
expansion timescale (or equivalently, the SNR age) depends upon 
the background ISM density, where low ISM densities facilitate a 
more rapid SNR expansion and vice-versa. For the ISM density in 
the range 0.1 to 1 particle(s) cm -3 , we find a corresponding SNR 
age of approximately 20 to 40 kyrs respectively. 

(ii) The magenta lines in Fig.[2]indicate the edge of SS 433's jet 
precession cone in the plane of the sky. By extrapolating the tra- 
jectory followed by the jet when in the plane of the sky, one would 
expect the jets to penetrate the W 50 nebula at the points P A , Pb> Pc> 
and P D (henceforth "the jet penetration points"), and alter the mor- 
phology of the nebula significantly. However, the radio observation 
of |Dubner et al.| ( |1998| > (see top panel of Figure |B3| > shows a sur- 
prisingly gradual morphological transition from the central circular 
region to the east and west lobes of the nebula, with no obvious ex- 
trusions at the jet penetration points. The lack of extrusions at the 
jet penetration points can be explained in one of to ways, either: 

[Scenario 1] - SS 433's jets have been somehow refocused back 
towards the precession axis, thus preventing any morphological de- 
formation at the jet penetration points. 

[Scenario 2] - SS 433's current jet precession state has not been 
active long enough for jets travelling at the current speed of 0.26c 
to reach the jet penetration points. 

(iii) The jet refocusing hypothesis has been investigated by var- 
ious authors (see for example Eichler (1983); Kochane k & Haw-| 
ley (1990); Peter & Eichler (1993)), and is a particularly appealing 
scenario because it could simultaneously explain the confinement 
of the X-ray emission (Brink mann et al.|2007| > to the region close 
to the jet precession axis. However, our hydrodynamical simula- 
tions ( Goodall et al. 2011 ) show that jet refocusing does indeed 
occur but via a different mechanism to those previously suggested. 
We find that refocusing of the jet trajectory occurs via momentum 
exchange with the turbulent interstellar gas contained within the 
nebula, which has been agitated by the jet itself. Although this refo- 
cusing mechanism does indeed provide an explanation for the con- 
finement of X-ray emission near the jet precession axis (see Figure 
8f and h of the companion paper (Goo dalTet al.|20lT] )), the side- 
effect of imparting momentum to the interstellar gas is to inflate the 
nebula, causing significant deformation at the jet penetration points 
(see Figure 8(e) and (g) of the companion paper ( [Goodall et al.| 
|2011| ). Thus, the refocusing of SS 433's hollow conical jets is un- 
likely to prevent nebular deformation at the jet penetration points 
marked on Figure|2] and we can safely rule out scenario 1. 

(iv) According to Figure [2] the nearest jet penetration point to 
SS 433 is P A at a distance of r A & 50pc. Assuming a jet speed of 
0.26c, the "in vacuo" jet travel time from SS433 to P A is ~600 
years. The jets do not travel in a vacuum however, and our hydro- 
dynamic simulations ( Goodall et al. 201 1 ) test the effects of various 
environmental properties (e.g. ISM density) upon the dynamics of 
jets with the same characteristics as SS 433's jets (precession angle, 
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jet speed etc). For comparison, simulating SS433's jets propagat- 
ing through an ISM density of 0.2 particles cm -3 , begins to have 
a noticeable impact on the SNR morphology at the jet penetration 
points approximately 3500 years after the jets are switched on. Dy- 
namical arguments can be used to estimate the likely age/duration 
of the current jet precession state at <3500 year^] This is consider- 
ably less time than the SNR age given in (i), and therefore favours 
the possibility of several jet episodes previous to the currently ob- 
served one. 

(v) Our hydrodynamic simulations (Go odall et al.||20l"T) also 
tested the effects of varying the jet precession cone angle. The re- 
sults indicate that any non-zero precession cone angles produce 
wider jet lobes than could explain the W50. Cylindrical jets are 
a special case of conical jets with a precession angle of zero, and 
only in this special case are lobes produced of the correct width 
when compared with the lobes of W50. Since SS433's current jet 
activity features precession with a cone angle of approximately 20°, 
it cannot be responsible for the lobes of W 50, and this is further ev- 
idence in favour of a previous jet episode having precession cone 
angle close to zero. 

(vi) Filaments A and J from Figure|2]lie directly along SS 433 's 
jet precession axis. The lack of movement from these filaments in- 
dicates an upper limit to the filament speed equal to the velocity res- 
olution of this investigation, of v res = 0.0405c. This speed is much 
lower than SS 433 's current jet speed of 0.26c, and the velocity dif- 
ference of 0.22c places a lower limit on the kinetic energy loss of 
jets, such that 85% of the kinetic energy per unit mass is transferred 
from the jets to the surrounding ISM, between jet launch and reach- 
ing the extent of filaments A and J. Figure[3]shows the velocity map 
from one of our hydrodynamic simulations of a persistently active 
jet with speed 0.26c. The velocity map shows that the persistent jets 
maintain very close to their launch speed even once the reach the 
full extent of W 50's lobes. The simulation featured an ISM density 
of 0.2 particles cm -3 which is a reasonable estimate for the ISM at 
the location of W 50 within the Galaxy. In fact, an untenably high 
ISM density would be required in the simulation, in order to ex- 
plain the reduction in jet speed to below v res = 0.0405c in the jet 
lobes. It is therefore unlikely that a persistent jet with speed 0.26c 
is responsible for the lobes of W 50. We propose two alternatives 
which would be consistent with the kinematics data: 

(a) A previous jet episode with a lower jet speed (i.e < v res ) could 
have created the jet lobes of W 50. This does not require an unfea- 
sible level of deceleration in the jets. 

(b) A previous jet episode of much shorter duration than the time 
required to produce the lobes, could be responsible. Upon cessation 
of the past jet activity, the transfer of kinetic energy from the jets 
to the surrounding ISM would then (and only then, without further 
input of energy from the power source SS 433) significantly decel- 
erate the jet lobes, in agreement with the data presented here and 
from our hydrodynamic simulations ( Goodall et al.|201l] ). 

Both of (a) and (b) require an additional jet episode dissimilar and 
previous to SS 433 's current jet activity. 

(vii) Intermittent jet activity is thought to be a common prop- 
erty of both microquasars and quasars ( Nipoti et al. 2005 ), and jet 
flares are routinely observed in other microquasars such as Cygnus 

6 Note that at higher ISM densities than 0.2 particles cm -3 the jets would 
require more than 3500 years to deform the nebula at the penetration points, 
and less than 3500 years at lower ISM densities, with a lower limit at the 
"in vacuo" timescale. 



X-3 (Mioduszewski et al. 2001 ; Miller -Jones et al.||2004] ), also 
GRS1915-105 ( [Mirabel & Rodriguez |fT994] > and V4641 ( |Hjellm-| 
ing 2000). For these microquasars however, the intermittency 
timescales are sub-year up to a few years, whereas we find hints 
of sub-millennium intermittency timescales for the SS 433 system. 



5 CONCLUSIONS 

In contrast to SS 433 's reputation as a persistent jet source, we pro- 
pose that the wealth of observational data on SS 433 and W 50 is 
consistent with SS 433 having experienced (at least) two distinctly 
different states of jet activity. First, the well- studied current jet state 
features t )a = 163 day precession with cone angle 0j et = 20.92° 
and jet speed Vj et = 0.2647c (Eikenberry et al. 2001 ) and is ob- 
servable both through the Doppler shifting of spectral emission 
lines and high-resolution radio interferometry. Second, a further 
state of jet activity which occurred some time in the past with the 
same axis of symmetry as the current jet system and either a zero 
(or very small) jet precession cone angle, or a precession period 
significantly longer than the duration of the jet activity, observ- 
able indirectly through the formation of the lobes of W50. We 
note that the observations of Cyg X-3's jet-flaring event in 1997 
were reported to have vi 991 > 0.81c and precession cone angle 
t/^1997 < 12° (Mioduszewski et al. 2001), and just four years later 
a speed v 2 ooi = 0.63c and precession cone angle i/^ooi = 2.4° were 
observed by |Miller J ones et al. (2004 ). The kinematic results pre- 
sented herein, underline the need for further investigation via hy- 
drodynamic simulations, of the possible intermittent jet activity in 
SS433's past. 
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Figure 3. We present a velocity map from high-resolution hydrodynamic simulations of the SNR-Jet interaction in SS 433 (see Goodall et al. ( 201^ for a 
detailed account of these simulations). A contour map of the Dubner-96 observation (white dashed line) has been overlaid upon the image to trace out the 
outline of the W 50 nebula. The above model can be summarised as a persistent cylindrical (zero precession cone angle) jet with speed 0.26c. The jet is incident 
upon an evolved SNR (blue circular region) of radius >40pc and age 21,000 years. The background density profile (no = 0.2 particles cm -3 , z = 40pc) has 
been optimised to match the observations of the east-west asymmetry in W50, and it takes 1,600 years for the jets to reach the extent of W50's lobes. The 
colour bar indicated the speed of the moving gas as a percentage of the speed of light, and the blue contour (solid line) indicates jet material which is travelling 
faster than the velocity resolution calculated in section [23] The harpoon-like light-pink/white regions (particularly evident in the eastern lobe) shows material 
which is travelling very close to the speed of the jet (0.26c) at the launch-point, which is a consequence of continuous ejection of kinetic energy from the 
persistent jet. 
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APPENDIX A: ADDITIONAL DETAILS REGARDING 
THE IMAGE PREPARATION 

Al Image Product Maximisation 

This method compares the positions of the 15 points sources in 
each image to accurately align the images. This works because the 
product of two similar 2D Gaussians is greatest when their peaks 
are coincident. 

A mask is first applied to each image, to set the image pixels to 
zero everywhere except for a small circular region around each of 
the point sources. Image product maximisation is then performed 
to align any two of the masked-images. 

The second masked-image I M 2 is shifted (resampled) by some 
small value dx in the x-direction, and dy in the y-direction. The 
images are then multiplied together and the sum of the product is 
recorder, according to: 

S (dx, dy) = ^ 7 M i (*, y) Imi(x + dx,y + dy) (Al) 

x,y 

The function S (dx, dy) is a 2D Gaussian (because the point 
sources being cross-correlated are also Gaussians). The (dx p ,dy p ) 
coordinate corresponding to the peak of S (dx, dy) denotes the shift 
that must be applied to the second image in order for both images 
to be optimally aligned. 



APPENDIX B: FITTING TO THE FILAMENTARY 
STRUCTURE OF THE W 50 NEBULA. 

The following images show the fits to the centroid positions of the 
filamentary structures common to each of the three VLA observa- 
tions. 
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Filament positions for the NVSS Observation. 
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Figure B2. The centroid positions are shown for the filaments of the NVSS-93 observation. 



© 0000 RAS, MNRAS 000, 000-000 



Kinematics of the interaction between the SS433 jets and the W50 nebula 




© 0000 RAS, MNRAS 000, 000-000 



